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SUMMARY 


The differential equations and boundary conditions appropriate to 
the analysis of torsional vibrations of hollow thin-walled cylindrical 
beams are presented. General solutions for the modes and frequencies of 
cantilever and free-free cylindrical beams of arbitrary doubly symmetrical 
cross section with uniform wall thickness are derived. Numerical calcu- 
lations for beams of rectangular cross section are used to determine the 
Influence on torsional vibrations of bending stresses due to torsion and 
of longitudinal inertia. The inclusion of longitudinal inertia exposes 
a strong coupling of the torsional vibration with longitudinal vibrations 
at low plan-form aspect ratios, and some of the frequencies obtained are 
seen to be actually associated with modes of primarily longitudinal vibra- 
tion. The solution for cylinders of rectangular cross section is also 
used to investigate the accuracy of a solution obtained from an analysis 
of an equivalent four-flange box beam. 


INTRODUCTION 


Investigations of the transverse vibration of hollow thin-walled 
cylindrical beams (see, for example, refs. 1 and 2} show that secondary 
effects such as transverse shear, shear lag, and longitudinal Inertia 
can have appreciable Influence on the transverse modes and frequencies 
of structures. Torsional oscillations are influenced by similar secondary 
effects, namely, longitudinal stresses which resist the tendency of the 
cross sections to warp out of their initial planes and the longitudinal 
inertia which Is developed by the warping motion. (The former effect Is 
commonly referred to as "bending stresses due to torsion" and will be so 
designated throughout this report.) 

In order to assess the importance of these secondary effects, a pro- 
cedure similar to that used in reference 1 has been used in the present 
paper to analyze the torsional vibrations of thin-walled cylindrical 
beams. The governing variational principle, differential equations, and 
boundary conditions are discussed, and general solutions for cantilever 
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and free-free cylinders of arbitrary doubly symmetrical cross section 
with uniform wall thickness are presented. ‘ Results of numerical calcu- 
lations based on the solution for the symmetrically vibrating free -free 
cylinder of rectangular cross section are used to illustrate the influ- 
ence of bending stresses due to torsion and longitudinal inertia on 
torsional oscillations . Finally, the results obtained from the calcu- 
lations for cylinders of rectangular cross section are compared with 
those obtained from an analysis of an equivalent four -flange box beam. 

SYMBOLS 

Ay cross-sectional area of flange 

A 0 cross-sectional area enclosed by median line of wall thickness 

Bj_ parameter defined by equation (29) 

C constant 

E Young's modulus of elasticity 

G shear modulus of elasticity 

Ip mass polar moment of inertia per unit length 

4A 

J torsional stiffness constant, —2 — 

7 T» 

K parameter, 

Xj 

Fourier coefficient 

L length of cantilever beam) half-length of free-free beam 

parameter defined by equation (36) 

T maximum kinetic energy 

U maximum strain energy 

a half -depth of rectangular beam 

b half -width of rectangular beam 
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jb m Fourier series coefficients 


i,j,m,n 


krji 


integers 

coefficient of longitudinal inertia, 


frequency coefficient, 



e./st 

l V i p 


p 

s 


t 

t' 


u(x,s) 

U F (x) 

X 



perimeter of cross section 

distance along perimeter of cross section (see fig. 1) 
thickness of effective shear-carrying area 
thickness of effective normal-stress-carrying area 
actual wall thickness 
displacement in x-direction 
displacement of flange in x-direction 
longitudinal coordinate 
shear strain 

longitudinal strain 


9(x) 

A 


P 

I 

P 


CT 

T 

0) 


rotation of cross section 
Lagrangian multiplier 
mass density of beam 

nondimensional longitudinal coordinate, ^ 

Xj 

distance from centroid of cross section to tangent to the 
median line of wall thickness (see fig. 1) 

longitudinal direct stress 

shear stress 

natural circular frequency of beam 
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BASIC EQUATIONS 
Assumptions 


The problem considered herein is that of the natural torsional 
vibration of a thin-walled hollow cylindrical beam. The cross section 
of the beam is assumed to be doubly symmetrical so that no transverse 
vibrations are introduced. The wall thickness of the cylindrical beam 
is assumed to be made up of an effective normal-stress-carrying thick- 
ness t', an effective shear-carrying thickness t, and a thickness t a 
associated with the mass properties of the beam. (These three distinct 
thicknesses are introduced for generality; they allow, for example, the 
inclusion of stringers which carry only normal stress and not shear.) 

For the present analysis, the following simplifying assumptions 
(also used in ref. 1) are made: 

(a) Changes in the size and shape of the cross section are 
negligible . 

(b) Stress and strain are assumed to be uniform across the wall 
thickness . 

(c) The effect of circumferential stress upon longitudinal strain 
is negligible . 

Because cross sections of hollow thin-walled cylinders are more 
easily sheared out of shape under torsional loads than under bending 
loads, assumption (a) has more significance in the present torsion prob- 
lem than it did in the bending problem In reference 1. For torsional 
vibrations, this assumption restricts the applicability of the analysis 
to structures thatmontain closely spaced bulkheads which are rigid in 
their own planes but offer negligible resistance to warping out of their 
planes. 

In agreement with the assumptions, the distortions of the vibrating 
beam will be defined by a rotation 0(x) of the cross section and a 
longitudinal displacement u(x,s) of each point of the median line of 
the beam wall. In agreement with the assumption of double symmetry of 
the cross section, the displacement u must be antiperiodic in the 
coordinate s with a period equal to one-quarter the perimeter p. 


Stress-Strain Relations 

The longitudinal and shear strains In terms of the displace- 
ments 9(x) and u(x,s) are 
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Sx 


and 


7 = ^ + 


P M 
M dx 


( 1 ) 


( 2 ) 


where p is the distance from the centroid of the cross section to the 
tangent to the median line (see fig. 1). The corresponding components 
of stress are 



and 


Txs " G (aS + 



( 3 ) 

(*) 


In elementary torsion theory, the cross section of the beam is con- 
sidered to be unrestrained against warping (that is, E = 0) with the 
result that the longitudinal stress 0 X is zero. 


Variational Principle and Energy Relations 
In order to use the variational equation 


5(U - T) = 0 (5) 

in the present analysis, the maximum strain energy U and the maximum 
kinetic energy T are expressed In terms of the displacements 0(x) 
and u(x,s). The variation Is taken Independently with respect to 0(x) 
and u(x,s) with the provision that these displacements satisfy the 
appropriate geometrical boundary conditions and that u(x,s) be anti- 
periodic in s with a period p/if. 

For a beam vibrating in a natural torsion mode, the maximum strain 
energy is 



The first term in equation (6) represents the strain energy usually con- 
sidered in elementary theory and the second term represents the strain 
energy contributed by the restraint of warping. 
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The maximum kinetic energy is 


T 



1 

2 



Cu 


£pt u^ds dx 


(T) 


where ca is the natural circular frequency of the particular mode, 1^ 

is the mass polar moment of inertia per unit length of the cross section, 
and p is the mass density of the "beam. The first term in equation (7) 
represents the kinetic energy usually considered in elementary theory; 
the second term represents the kinetic energy of warping or longitudinal 
mot ion . 


Differential Equations and Boundary Conditions 

The substitution of equations (6) and ( 7 ) into equation ( 5 ) and the 
application of the calculus of variations results in the simultaneous 
integrodifferential equations for the displacements u(x, s) and 0(x), 


Et* 


6^u + q 6t chi 
6x 2 6s 6s 


+ Gt 


6^u + q d(pt ) as 

3 s 2 6 s dx 


i-LOD^tgU 


0 


(8) 



+ P 



ds + IpO) 8 = 0 


( 9 ) 


together with the boundary conditions at each end of the beam. 


and 



( 10 ) 


(H) 


At a fixed end, both boundary equations (10) and (ll) are satisfied 
since the geometrical boundary conditions require that both Su and 58 
be zero.' At a free end, since both 6u and 58 are arbitrary, equa^" 
tions (10) and (ll) yield the natural boundary conditions 
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and 



(13) 


Equation (12) expresses the condition of zero longitudinal direct stress 
and equation ( 13 ) expresses the condition of zero net torque on the free 
end. 


Although the integr ©differential equations (8) and ( 9 ) were obtained 
by use of the variational principle, these equations can, of course, be 
obtained directly from consideration of equilibrium conditions. 


GENERAL SOLUTIONS FOR CYLINDERS OF UNIFORM WALL THICKNESS 


In this section, the exact solutions are presented for cylinders of 
uniform wall thickness; that is, t’, t, and t a are each uniform around 
the cross section. These solutions are obtained by means of Fourier 
series and the variational principle (eq. ( 5 )) - a more expedient approach 
than a direct attack on the integrodifferential equations ( 8 ) and ( 9 ) and 
their associated boundary conditions. 


Cantilever Beam 

For a cantilever beam (fig. 1(c)), the geometrical conditions are 
u(0,s) = 0(0) = 0. The natural boundary conditions, equations (12) 
and (13) , are satisfied as a result of the variational process. The 
displacements 9(x) and u(x,s) are assumed to be given by the 
expressions 


and 


0(x) = C + 


m=l, 3 > 5 


, mrtx 

b_ cos 

^ 2L 


u 


mrtx 2nrts 


(x,s) = y a sin sin 

' ’ ’ ■!— , mn ot 


m=l,3,5 n=2 ,4,6 


2L 


( 14 ) 


(15) 


The choice of the particular trigonometric functions used in the expres- 
sions for u and 0 was guided by consideration of the orthogonality 
required for the simplification of the strain-energy expressions. The 
condition u(0,s) =0 is satisfied by each term of equation ( 15 ). The 
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constant C in the expression for 0(x) is required to allow 9(L) to 
he unrestricted. The condition 


0(0) = C + 21 b m = 0 (16) 

m=l,3,5 

is introduced into the variational procedure by means of the Lagrangian 
multiplier method. 

Substitution of equations (l4) and (15) into equations ( 6 ) and ( 7 ) 
yields 



In order to make equation (17) stationary and at the same time 
satisfy the constraining relationship 
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cp = C+ 2_ ^ = 0 

m=l,3,5 


(18) 


it is sufficient to set 


S(U - T + Acp) = 0 


(19) 


where the variation is with respect to the a's, b’s, and C considered 
as independent variables and where A is a Lagrangian multiplier. The 
variational process yields the following equations: 


$ 2 - 1*1 = ft.tef a a , + Gt (&L) £B 

Baj^j \2L/ 4 \ p / 4 


^a T a U = ° 


- Gt *f- K J b i - 


i - 1, 3, 5, . . 

j = 2 , k, 6 , . . 


( 20 ) 


3b! &>i k ^ (slJ k n= fj- j6 "■> 1 


i -1 


Gt (lf f ^ ^ i bi ' (_1) S “V + * 


(i = 1, 3, 5. ■ • •) (21) 


+ A 

ac ac 


^ mrf.,3,5 ^ 


m -1 

g(-l) b,,, - co 2 IpLC + A 


= 0 ( 22 ) 


where 



10 


NACA TN 3206 


Kn = |^P cos ^p. ds (23) 

and 

“o - as - ^ (24) 

in which A 0 is the area enclosed by the median line of the wall thick- 
ness. In obtaining equation (21) the relation 

i <£p 2 ds = Ko 2 + | XI K^ 2 ( 25 ) 

y J n=2,h,6 


was utilized. This results from an application of Parseval's theorem, 
Kq and being coefficients of a Fourier series expansion of p. 

With the use of the nondimens ional parameters 




!pL 2 2 


GJ 


cu 


(2 6) 




(27) 


K 2 ^ 

GL 2 


(28) 


and 



( 29 ) 


equations (20), (21), and (22) may be reduced to 
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(Bi 2 + W -0 I 1 ' 1 ’ l’ l ) (30) 

\j = 2, K, 6, . . ./ 


Lt ? ^ — 

* * 6 ^ 


ain+ ^-i 1+ *nS-,6 K “ 2 )-i kr2 ' ) i 


2 1 ^2 


— (-1) bj,^ + — = 0 

irt 1 GJ 


(1 = 1, 3, 5, . . •) (31) 


2_ ^ (-1) 2 Uta + % 2 C - g = 0 

m=l,3j5 x GJ 


where J is the elementary torsional stiffness constant defined as 


i^Ao^t 


Fran equation (30), 


^ L K j lJ 

1J = Bi 2 7 16 J 2 "i 


i - lj 3> '?) • • • 

<5 = 2, 6, • • • 


Substituting this expression for a. i into equation (31) and solving 
for bjL yields 


hi = 


— (-1) 2 kJc - — 

ijc k ; TC L GJ 


(i = 1, 3, 5, • • •) (35) 


ptB i £ 


n=2 ,b,6 B i + l6n 2 / 2 


.lv 2 

2 x 
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The substitution of equation (35) into equation ( 32 ) and into the 
constraining relationship (l 8 ) yields the following two homogeneous 
equations in C and A: 


V 


1 + 


^ (=) i 

m=l,3,5 ™ 


C - 


1 + 




m=l,3,5 


_2_ 

mrt 


m-l" 


(-D^ 

AL 

N_ 

"m 

GJ 


£ 9=0 (37) 


00 ~ 

1 4. k 2 y~ — 

m-1 

(-D 2 

C - 

( II 


1 + %! 4^ _ mrt 

m=l,3,5 



l m=l,3,5 

N m J 


(58) 


The condition for a nontrivial solution for C and A gives the 
frequency equation 


*T 


1 + krji 


m=l,3 >5 


(~) " 


1 + 


V 


m=l,2,3 


_ 2 _ (- 1 ) 

urn 


m -1 

2 


0 2 (-1) 


m-1 

2 


1 + ^T ^-= _ mit N 

m=l,3>5 m 


CO 

m=l73,5 % 


= 0 (39) 


which the frequency parameter hp must satisfy. Each tern of each of 
the- infinite series in-the frequency equation contains therefore, 
the roots of equation ( 39 ) are most readily obtained by trial. 


Once kip has been determined, a 
by letting C = 1 and solving either 
then finally evaluating b . and a., ± 
and (34). J 


particular mode shape can be found 
equation (37) or ( 38 ) for A, and 
successively from equations (35) 


Free-Free Beams 

Symmetrical modes .- For a free-free beam of length 2L with the 
origin at the midspan (see fig. 1), the form of the Fourier series 
assumed for 0 (x) and u(x,s) when the beam is undergoing a symmetrical 
mode of torsional vibration may be the same as that assumed for the 
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cantilever beam, of length L. (See eqs. (l4) and ( 15 ).) However, if 
consideration is given to the right half of the beam, the only geometri- 
cal boundary condition enforced is u(0,s) = Oj therefore, the con- 
straining condition (l8) does not apply and the appropriate frequency 
equation is obtained from equation (37) by letting X = 0. The frequency 
equation is 



(40) 


After a particular root top has been found from equation (4o), the 

corresponding mode shape may be obtained from equations (35) and (34) 

(with A = 0) . 

Antisymmetrical modes .- For a free-free beam of length 2L vibrating 
in an antisymmetrical torsional mode, specific consideration need again 
be given only to the right half of the beam (see fig. 1). With the ori- 
gin taken at the beam midspan, the only geometrical boundary condition 
imposed on the right half -beam is that 0(0) = 0. The longitudinal dis- 
placement u(x,s) is unrestricted over the length of the beam. 

The displacements 0(x) and u(x,s) for the antisymmetrical modes 
are assumed to be given by the expressions 


and 


8(x) 


GO 

Cx + 5 

m=2,4,6 


bjn sin 


mrtx 

2L 


u (x, s ) 


m=0,2,4 



Smn cos 


mjtx 

2L 


sin 


2mts 

P 


0U) 


(* 2 ) 


As before, the choice of the particular trigonometric series in 
equations (4l) and (42) was guided by considerations of orthogonality. 
The linear term Cx appearing in equation (4l) is necessary to allow 
sufficient freedom for the beam to twist. 


Using equations (4l) and (42) in equations (6) and ( 7 ) gives 
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U - T. 


0 J \m=2,4,6 n-2 ,4,6 


mrt . mrtx . 2mts 

— sin — — sin ds dx + 

2L ^ nn . 2L p i 


Cf 


m=0,2,4 n=2 ,4,6 


2nrt _ rmtx 2ruts , 

-p- aim cos 2T COS —p~ + 


CO 

mjt . mrtx , _ 

45,6 ^ cos ^ + pc 


ds dx - 


ipio 2 r L 


/ 0X + £ ^ 

J 0 \ m=2,4,6 


. mrtx \ , 

sin dx - 


p.t a cu‘ 


>2 P L 


0 J \m=0,2,4 n=274,6 


mtx . 2nits . , /i, _v 

a nin cos sin ds dx (43) 


The variation of equation (43) vith respect to the a's, Vs, 
and C gives, after simplification. 


( B i 2 + + k | iJKjbi = 0 ^ 


i = 2, 4, 6, . . . 
j “ 2, 4, 6, . . . 


g * 2 ^ inVin + ^ P -( 1 + Zy ^ ^ 

^ n=2 ,4,6 |_ 8 \ 2,1 n=2,4,6 


2 ) - | *T 2 *i - 


4- (- 1 ) 2 % 2 LC = 0 
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From equation (44), 



±iSL b 

Bj 2 + 16 3 2 


/i - 2 , 4, 6 , . . A 
\J = 2 , 4, 6 , . . ./ 


(48) 


which, except for sign, is the same expression as that obtained for the 
cantilever and symmetrically vibrating free-free beams (eq. (34)). From 
equation (46) , 


a 0j “ 


8 P K . 
* L V 


LC 


(j = 2 , 4, 6 , . . .) (49) 


gives 


16 j 2 + B q 2 

Substituting equation (48) into equation (43) and solving for b^ 


I ^ 

iN, 


LC 


(i = 2 , 4, 6 , . . .) ( 30 ) 


where is given by equation ( 36 ) . 
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Substitution of equations (49) and (50) into equation (47) and 
simplification gives the following frequency equation for the antisym- 
metrical vibrating free-free beam: 



(51) 


After a particular value of kj is obtained from equation (51) , "the 
corresponding mode shape may be computed by letting C = 1 and calcu- 
lating the b's and a's successively from equations (50), (49), 
and (48) . 


Discussion of Parameters 

In the frequency equation, the unknown natural frequency appears 
only in the frequency coefficient kp which is defined by the commonly 

used relationship a) = 

the effect of bending 
be neglected, it is sufficient to set K equal to zero in the final 
frequency equation. The parameter k TT is associated with the effects 

of longitudinal Inertia and appears only in the expression for B i . 

(See eq. (29)-) If the effects of longitudinal inertia are to be neg- 
lected, k LI may be set equal to zero. With k^. = 0, B_ L is independ- 
ent of kr£ and simplification of the trial-and-error solution for the 
natural frequencies results, since the infinite summation contained In 
% (eq. (56)) need be calculated only once for a particular beam. In 

the following section, it is shown that the influence of longitudinal 
inertia is often negligible. 

The K^'s are the Fourier coefficients of the function p(s), which 
depends only on the contour of the cross section of the beam. It Is 
interesting to note that, for a cylinder whose cross section is either 
circular or a regular polygon, the K^'s become zero and the frequency 
equations (59), (40), and (51) can be written In a closed form identical 
with the frequency equations obtained from elementary theory. ThiB Is 
consistent with the generally known fact that beams of such cross section 
have no tendency to warp; hence there can be no longitudinal inertia or 
bending stresses due to torsion. 


. The parameter K is associated with 


stresses due to torsion, and. If this effect is to 
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RESUITS FOR CYLINDRICAL BEAMS OF RECTANGULAR CROSS SECTION 


In order to illustrate the effects of bending stresses due to 
torsion and longitudinal inertia on the natural torsional frequencies, 
the present analysis has been applied to hollow thin-walled cylindrical 
beams of rectangular cross section. The be ams considered are assumed 
to have walls of constant thickness, with t* = t = t a , and a cross- 
sectional aspect ratio of 3*6. The numerical calculations have been 
limited to the free-free symmetrical modes of vibration. A value of 
e/g of 2.65 (appropriate for aluminum alloys) was used. Results in the 
form of the variation of the frequency coefficient kp with plan-form 

aspect ratio L/b are shown in figure 2 for the first four modes of 
vibration. 

The short-dashed curves in figure 2 show results obtained from 
calculations based on elementary torsion theory where the restraint due 
to warping of the cross section Is neglected. These curves serve as a 
basis for resolving the effects of bending stresses due to torsion. The 
long-dashed curves are calculated from equation (40) with k LI = 0 and, 

consequently, represent the natural torsional frequencies when only the 
effects of bending stresses due to torsion are included. (For the rec- 
tangular tube the quantity N-p that appears In equation (40) can be put 
in the closed form presented In appendix A; this closed form was used in 
this calculation as well as In the calculation of the exact frequencies 
mentioned In the next paragraph.) Comparison of these two sets of curves 
shows that , as might be expected, the effect of including bending 
stresses due to torsion is to increase the calculated natural frequency 
of the beam. For the higher plan-form aspect ratios, the contribution 
of bending stresses due to torsion is small. However, this contribution 
does become appreciable at the lower values of L/b for the higher modes 
of vibration. 

The solid curves show the calculated results for the first four 
modes, when the combined Inf luence' of longitudinal inertia and bending 
stresses due to torsion is taken Into account. The values of kp shown 
by these curves, which are referred to as "exact," are the lowest four 
roots of equation (40). 

Comparison of these curves with the long-dashed curves shows that, 
as should be expected, the effect of including longitudinal inertia is 
to decrease the calculated torsional frequencies. For plan-form aspect 
ratios L/b above 4, this effect Is practically negligible. For the 
lower plan-form aspect ratios, however, the effect becomes appreciable; 
the exact curves no longer follow the corresponding long-dashed curves, 
but drop off in a seemingly disorderly manner. 



18 


NACA TN 3206 


The explanation for the apparently erratic behavior of the exact 
curves may be obtained by considering the kinetic” energies included in 
the derivation of the exact frequency equation (eq. (LO)). As was pre- 
viously mentioned, the kinetic energy (eq. ( 7 )) consists of two parts: 
one which represents the contribution of the rotational inertia and 
another which represents the contribution of the longitudinal inertia. 

If only the inertia of rotation is considered, the resulting frequency 
equation represents a primarily torsional vibration in which only the 
effects of bending stresses due to torsion are included. (The results 
of such a solution have already been shown by the long-dashed curves in 
fig. 2.) If, however, only the longitudinal inertia is considered, the 
frequency equation obtained represents a particular type of primarily 
longitudinal vibration which is hereinafter referred to as "warping" 
vibration, (it should be noted that this vibration is not entirely 
longitudinal; some twist is present, just as longitudinal motion is 
present in the torsional vibration.) The exact frequency equation, 
which includes the effects of both the rotational and the longitudinal 
inertia, can be looked upon as representing the coupling of the two 
"uncoupled" vibrations (torsional and warping). 

In order to illustrate this coupling, the frequencies of the warping 
vibration have been obtained In appendix B and are included in figure 2 
as the dot-dash lines. Consider now, for example, the fourth exact fre- 
quency. As the plan-form aspect ratio is_decreased, the exact frequency 
follows the fourth uncoupled torsional frequency until an L/b value of 
about 3-5 is reached. At this point, the effect of the warping vibra- 
tions becomes predominant and the curve falls off rapidly until the third 
uncoupled torsional frequency Is reached. (Presumably, the mode shape of 
the vibration during this transit between the uncoupled torsional fre- 
quencies is roughly the same as the first warping mode.) As L/b is 
decreased further, the curve continues to cascade in a similar manner. 

Thus the behavior of the exact frequency curves is seen to be quite 
rational. At a given value of L/b, the exact frequency coefficients 
and their associated modes can be interpreted by considering the fre- 
quency coefficients of both uncoupled torsional and warping vibrations. 

At L/b = 2, for example, the first, second, and fourth exact frequencies 
are comparable, respectively, to the first, second, and third uncoupled 
torsional frequencies. The third exact frequency, which appears to be 
an extra torsional frequency, Is associated with the first warping mode 
of vibration. This same phenomenon of seemingly extra frequencies due 
to inclusion of -longitudinal Inertia was noted for bending vibrations by 
Traill -Hash and Collar in reference 2. 

The disposition of the exact curves in figure 2 with respect to the 
approximate curves suggests that it may not be necessary to consider the 
problem of coupled torsion and warping in order to obtain natural fre- 
quencies of reasonable accuracy. Certainly, for the higher values of 
L/b, the frequency coefficients obtained by neglecting longitudinal 
Inertia are sufficiently accurate. For the lower values of L/b, the 
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transitions of the exact frequencies can be characterized by considering, 
in addition, the frequencies of uncoupled warping vibration. 

Besides showing the quantitative influence of the secondary effects, 
the present analysis can be put to another practical use: to assess the 

accuracy of simplified torsional-vibration theories. One such theory is 
presented in appendix C and is based on an analysis of a four-flange box 
such as is shown in figure 3* In such a box beam, the flanges are 
assumed to carry only normal stresses while the sheets carry only shear. 

In order to represent the rectangular cylinder as a four-flange beam, 
the cross-sectional area of each flange is taken to be equal to l/6 the 
section area of the walls adjacent to each corner of the box. In the 
analysis, the effects of longitudinal inertia are neglected. 

The results of calculations for the natural torsional frequencies 
of rectan gular tubes based on the four-flange solution are shown in 
table 1. For comparison, the frequency coefficients calculated from 
the elementary solution and from the exact solution without longitudinal 
inertia are included. Examination of the results shows that the four- 
flange solution gives frequencies that are in good agreement with those 
of the exact solution, even for the low-aspect -ratio box beam where the 
influence of bending stresses due to torsion is important. 

CONCLUDING REMARKS 

The results of numerical calculations show that elementary theory is 
adequate for the analysis of torsional vibrations of box beams with high 
or medium plan- form aspect ratios. For box beams with low plan-form aspect 
ratios, the influence of bending stresses due to torsion becomes important. 
The Influence of longitudinal inertia, however, does not become important 
until extremely low (and probably impractical) plan-form aspect ratios are 
reached. A simplified procedure for obtaining torsional frequencies, based 
on an analysis of an equivalent four-flange box beam, was found to give 
adequate predictions of the effects of bending stresses due to torsion. 

The general analyses presented, as well as the numerical results for the 
rectangular box beams, should be useful in the assessment of the accuracy 
of other simplified procedures that may be developed. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va., May 12, 195^* 
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APPENDIX A 

CLOSED FORM OF N ± FOR CYLINDRICAL BEAMS 
OF RECTANGULAR CROSS SECTION 


The expression for N^ can he written in closed form for a cylin- 
drical beam of rectangular cross section and constant thickness. For 
such a beam with the dimensions shown in figure 2, the parameters become 


p = 4(a + b) 
j - l6a^b% 

a + b 


Kn 


4 a - b 
n n 


sin 


2mtb 


(n even) 


*n = ° 


(n odd) 


(Al) 


With equation (Al), equation ( 36 ) can be written as 
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The Infinite summations in equation (A3) 
form and the resulting closed expression for 


can be written in closed 
Nj, is 


N i = 


(i *)‘ 


1 + 


'a 2 - b 2 ^' 


ab 


ab 


\4(a + b)‘ 


1 


a - b 

cosh — rr cosh 

8 a + b 


JtB. 


_ 1 v 2 

2 TC 


(A4) 
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APPENDIX B 

SOHJTION FOR THE WARPING VIBRATION OF CYLINDRICAL BEAMS 


The solution for the warping vibration of thin-walled cylindrical 
beams can be obtained from the variational equation, as was mentioned 
previously, by setting the first term in the expression for kinetic 
energy (see eq. ( 7 )) equal to zero. 

In order to obtain the anti symmetrical warping modes that couple 
with the free-free symmetrical torsional modes of cylinders with uniform 
wall thickness, it is sufficient to set Ip equal to 0 in equation (17)? 
then 


* Jo J \m=l,3,5 n=2 ,4,6 ^ 


cos 3!52£ sin ^ 3ts A ds dx + 
2L Pi 


Gt 


oL „ 

^0 ^\p=l75,5 n=2,k } 


2mt _ mitx 2mts 

■ 1 sin cos 

g p 2L P 


p £ ^• b m sln lr) d6te - 

m=l73,5 211 2L / 




/W £ 

^0 ^ \m=l,3>5 n=2 ,4,6 


mrtx 2 m s\ 
p / 


2L 


(Bl) 


The substitution of equation (Bl) into equation ( 5 ) yields the following 
equations: 


(b ± 2 + 16 3 % - k 2 K 3 i 3bi - 0 (j : £ J; 2; ; : ;) 

i j 


(B2) 
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d(U - T) 
dbj 


— it 2 inK n a in + - 

n=2 ~k,6 K ^ ^ 8 


+ * fl < 

23 n=2 X6 


bi = 


(i = 1, 3, 5, • . 0 (B3) 

If equation (B2) is solved for and the resulting expression is 

substituted into equation (B3 ) , the following equation is obtained: 


( j ») 2 

8 



£ 

=9 li 



M y- n V V 

J n=2 ~£,6 B 2 + l6n 2 J 1 


= 0 


(BIO 


In order that nonzero values of 
equation 


bj_ may exist, the frequency 



K 


_EL 


,6 B^ + i6n £ 


= 0 


(B5) 


must be satisfied. The frequency coefficient foj is contained In the 
parameter B^_ . Once a particular root bp has been found from equa- 
tion (B5) , the corresponding mode shape can be obtained from equation (B2) 
by setting bq equal to the arbitrary value of 1. 
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APPENDIX C 

SOLUTION FOR A FOUR-FLANGE BOX BEAM 


A simplified torsion theory, used to obtain the torsional modes 
and frequencies of- structures where the influence of bending stresses 
due to torsion is to be included, can be based on an analysis of a four- 
flange box beam such as that shown in figure 3* (The torsional vibra- 
tions of such a structure were considered in 1951 by E. H. Mansfield of 
the Royal Aircraft Establishment, who obtained results for a cantilever 
beam and also for an antisymmetrically vibrating free-free beam; in this 
appendix results will be shown for the free-free symmetrical modes.) In 
an analysis of such a structure the following assumptions are made: 

(a) The flanges are assumed to carry only normal stresses. 

(b) The sheets connecting the flanges are assumed to carry only 
shear . 

(c) The cross sections remain undistorted. 

(d) Longitudinal inertia is assumed to be negligible. 

The actual flange of the box beam should be replaced by an effective 
flange which incorporates all of the normal-stress carrying capacity of 
the beam. 

In agreement with these assumptions, the distortion of the vibrating 
beam is defined by a rotation 9(x) of the cross section and a longi- 
tudinal displacement Up(x) of each flange. 

The longitudinal and shear strains are given by equations (l) and ( 2 ) 

with 

Su = 

Ss "b 

and 

p = a (C2) 


(Cl) 


for the cover sheets and 
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and 


p = b 


(C4) 


for the webs . 

The solution for the four-flange box beam could be obtained by 
means of the variational principle; however, the direct use of the 
integrodifferential equations is more expedient. 

Substituting equations (Cl) to (C4) into equation ( 9 ) yields 


-tot (b - a)up' + totab(a + b)e" + a> 2 I p e = 0 (G5) 

where primes denote differentiation with respect to x. Equation (C5) 
can be recognized as the equilibrium equation for a cross-sectional 
element of the beam. 

Because of the assumption that the sheet in a four-flange box 
carries only shear, equation (8) has meaning only at the flanges of the 
box. Since the flanges are assumed to be concentrated at a point, the 
thickness t' in equation (8) must be defined as 

t' = 8(s 

where 8 is the Dirac delta function 

8(s - b) 

8(s - b) 

in such a way that 



Therefore, if equation (8) is integrated over the corners, that is, from 
s=b-A to s = b + A, where the infinitesimal quantity A approaches 
zero, the following differential equation is obtained: 

EA-pUjp" - Gt^i + ^Up + Gt(b - a)0* =0 


- b) (C6) 

defined by 

= 0 
= 00 


(s ± b) 
(s = b) 


(C7) 


(C9) 
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Equation (C9) is the equilihrium equation for an element of a flange. 

In order to obtain the equilibrium equations (C5) and (C9), use was 
made of the previously derived integrodifferential equations ( 8 ) and ( 9 ) 
to show the applicability of these general equations to a particular 
problem. However, equations (C5 ) and (C9) could have been written 
directly from consideration of the equilibrium of the various elements. 

Either 0 or u^, may be eliminated from equations (C5) and (C9) 
to give the following equation in 0 or Up alone: 




4ab 

(a + b )2 


kGtL 2 


EA-o(a + b) 


d 2 0 

d| 2 


" V 


kGtL 2 


EAp(a + b) 


0 = 0 


(CIO) 


or 


d 4 up 

< 1| 4 




4ab 




(a + b ) 2 ®Ap(a + b) 


d % 


d| 


P 4GtL 2 

** —7 rr ^ 

EAp(a + b) v 


= 0 (Cll) 


where 

Solutions of equations (CIO) and (Cll) are of the form e a ^ where 
a is a root of the equation 


a 4 + top 2 —Jiab. 

(a + b 


4GtL 2 


b) 2 + "b)J 


Wtff* b) " ° (C12) 


The general solutions for 0 and up are given by 


0 = C-p sinh Ojt + Cg cosh + Cj sin ctgl + C^ cos Og| (C13) 


and 


Up = Dj_ sinh a-pi + I >2 cosh aq| + sin agg + D 4 cos agg (Cl4) 
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where cx^ and are the real and imaginary roots, respectively, of 

equation (C12) . The constants Dj_ to are related to the con- 

stants Cq_ to C 4 . through equation (C5) or (C9) as follows: 


where 


L 2 

b 2 - a 2 

4 

a 2 b 2 

L 2 

b 2 - a 2 


a^b^ 

L 2 

b 2 - a 2 

if 

a 2 b 2 

CM 

b 2 - a 2 

4 

a 2 b 2 


°1 

“l 

02 

02 



+ 


+ 


*I 2 ) C 2 

k>t 2 )c 1 

)ci f 

% £ ) c 5 



= 0 


► 


= 0 



(a + b ) 2 
ifab 


(C15) 


(C16) 


In addition to equations (C15), four other relations between the 
constants are needed and these relations are obtained from the boundary 
conditions imposed on the beam. 


For a free-free beam undergoing symmetrical vibrations, the boundary 
conditions are at the midspan 


*£ = 0 

dx 


” r_0 J 

and at the free end 


dup 

dx 


(C17) 


(Cl 8 a) 


(a - b)u„ + ab(a + b) — = 0 
* dx 

where equation (Cl 8 b) is the condition of zero torque. 


(Cl 8 b) 
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The conditions (C17) and (Cl8) can be expressed in terms of either 
set of constants^ by the use of equations (C15) . The condition for a 
nontrivial solution gives the frequency equation 


(Bag 2 - lSj 2 ^ 2 


^Ba ^ 2 + bp 2 ja 


cosh sinh oq_ ag cos ag -ag sin ag 

(Vj 2 - Beiges Inh aj_ - Ba^coah ^Sj 2 + sin ag ^Icj 2 + Bcq 2 ^ cos ag 


0 (C19) 


which reduces to 


(^T 2 + cosh sin ctg - - Ba^Jag sinh cos => 0 (G20) 


The frequency equation has been obtained for the special case of a 
rectangular box beam with web and cover sheets of equal thickness in 
order to compare the numerical results from the four-flange solution 
with the results from the exact solution. However , the solution for 
the four -flange box beam with unequal web and cover sheets or with other 
boundary conditions may be obtained by means of the same analysis 
procedure . 
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TABLE 1 


COMPARISON OF RESUITS OBTAINED FROM FOUR-FLANGE SOLUTION WITH 


THOSE OBTAINED FRCM ELEMENTARY AND EXACT SOLUTIONS 



L/b 

Elementary 

Exact solution without 

Four-flange 

solution 

longitudinal inertia 

solution 

First mode top 

2 

3.14 

3.30 

3-31 

6 

3.14 

3.18 

3-18 

10 

3.14 

3.16 

3.16 

14 

3.14 

3.15 

3.15 

Second mode top 

2 

6.28 

6 . 9 8 

7.11 

6 

6.28 

6.53 

6.54 

10 

6.28 

6.39 

6.4o 

14 

6.28 

6.35 

6.35 

Third mode 

2 

9.42 

— 

10.80 

11.01 

6 

9.42 

10.02 

10.13 

10 

9.42 

9- 78 

9-77 

14 

9.42 

9.65 

9.63 






1 


l 

y 


1 

y 


(c) Cantilever beam. 

(d) Free-free 

beam. 

9* 

Figure 1. — Coordinate 

systems and sign 

conventions. 
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